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Recent theoretical ideas and observational claims suggest that the fine structure constant a may be variable. We examine 
a spectrum of models in which a is a function of a scalar field. Specifically, we consider three scenarios: oscillating a, 
monotonic time variation of a, and time-independent a that is spatially varying. We examine the constraints imposed upon 
these theories by cosmological observations, particle detector experiments, and "fifth force" experiments. These constraints 
are very strong on models involving oscillation, but cannot compete with bounds from the Oklo sub-nuclear reactor on models 
with monotonic time-like variation of a. One particular model with spatial variation is consistent with all current experimental 
and observational measurements, including those from two seemingly conflicting measurements of the fine structure constant 
' using the many multiplet method on absorption lines. 
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The conventional view is that the fine structure constant a, the coefficient that determines the strength of the 
electromagnetic interaction, is a constant. There is an effective change in a at energies greater than 1 GeV due to 
renormalization effects associated with electroweak symmetry breaking. However, at temperatures far below the GeV 
scale, and certainly since matter-radiation equality, these effects are negligible. 

There are recent claims of a statistically significant variation of the fine structure constant at large redshifts 
> I 0.2 < z < 3.7 (cf. [1]), 

'" - 2< ' a < (8 3 f ) ;° ( " ^ (-0. 5 7 ± 0.10)x 1 0-', (1, 
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where these measurements are averaged over quasar absorption systems in the given redshift range. An independent 
measurement using similar techniques is reported in [2], 

^ I a(0.4<z<2.3)-a(z = 0) = (-0.06 ± 0.06) x 10-. (2) 

*j-H a(z = 0) 

One issue discussed in this paper is whether these two results are contradictory. The only other possible statistically 
significant detection of variation of a coupling constant comes from measurements of isotopic abundances from the 
sub- nuclear reactor at the Oklo mine, where the analysis allows two possibilities for the value of a [3-5]. The favored 
branch is a null result providing a strict bound on variation, but the other branch would also yield statistically 
significant variation. A full discussion of experimental bounds on the variation of coupling constants is given in [6] . 

There is also considerable theoretical motivation for the possibility of variable coupling constants. Coupling con- 
stants appear to be naturally variable in unified theories [7] , and among other things will vary with a changing dilaton 
[8]. A changing fine structure constant would be accompanied by variation in other fundamental parameters, including 
the grand unification scale [9] . Variation in the fine structure constant also can be associated with violations of local 
Lorentz invariance and CPT symmetry [10]. 

In this paper, we consider experimental and observational constraints on field theories allowing variation of the 
fine structure constant, as well as whether any simple models satisfy both these constraints and the measurements of 
Webb et al. In § III, we consider whether finite temperature field theory provides a sufficient correction to account for 
these observations. In the remainder of the paper, we consider three classes of models and discuss phenomcnological 
constraints that must be imposed. Models with predominantly oscillatory behavior in the recent past are discussed 
in § IV. Models with monotonic variation are considered in § V. Finally, spatially varying but time-invariant models 
are considered in § VI. For the remainder of this paper, we will use spatial variation to refer to models in which the 
value of a is different when measured at two points separated by space-like vectors and similarly time variation will 
refer to differences between points separated by time-like vectors. While not every model will fall strictly into one of 
these three categories, the goal here is to produce representative constraints applicable to a wide range of models. 
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II. GENERAL CONSIDERATIONS 



If a coupling constant is variable, we can treat it as a function of some field. This paper will consider models in 
which the fine structure constant is a function of a scalar field. Clearly not every possible scenario can be described 
in terms of a scalar field, but for small changes in a a scalar field will cither be exact or a good approximation for a 
very wide variety of models. Models that cannot be described in terms of a scalar field fall outside the scope of this 
paper. 

The QED Lagrangian including a real scalar field <j> is 

Cqed = \(d^) 2 - V(<j>) + l(i^d» + erf - m )l-± (0„(eA,) - d^eA^)) 2 , (3) 

where I is summed over all leptons I = e,fi, r, and A^ is the vector potential. Note that for constant e and under the 
rescaling eA^ — > A^, (3) becomes the more familiar 

Cqed = \{d^f - V{4>) + + 7M„ - m,) I - V"- ( 4 ) 

Both forms of the Lagrangian have been given simply to illustrate the proper transformation that leaves the two 
entirely equivalent. The form given in (3) is given as the primary form for ease of calculation. 
If a has ^-dependence, the QED Lagrangian will include a term 

Cqed = ■■■ + ly/^a^'fA^l. (5) 

The equation of motion for a scalar field <f> in an expanding universe is 

'<j> + 3H(t)4>=-V'{4>), (6) 

where V denotes the derivative with respect to <f>. 

As an illustrative example, consider the case in which cj> has a minimum at <j>o = 0, and the potential V(<f>) is the 
harmonic potential 

V^^^ml^-M 2 . (7) 

Substituting, the equation of motion becomes 

4> + 3H (t)<h + ml(<j> - <Ao) = 0. (8) 

This is the equation for a damped harmonic oscillator, with damping term 3H(t) (which has units of time -1 ), and 
therefore admits the following characteristic solutions today: 

1. Large m^. If rn^ » H 7 cj), and therefore a, will be rapidly oscillating with period T <C 3H 0} where H is the 
value of the Hubble parameter today. The universe out to a redshift of 3.7, where variation of the fine structure 
constant has been reported in Webb et al., is observed over a time period of approximately 1/Hq. Therefore, 
over our range of observations, these models would predict that a would appear to be oscillating. Models with 
these features are discussed in § IV. 

2. Intermediate m^. If ~ H, <fi, and therefore a, could be oscillating with period T ^> 3H . In this case, we 
should observe the variation of a as monotonic rather than oscillating over our range of observations. Also note 
that the damping term yields a characteristic decrease in a of one e-folding per 1/3H - Models with monotonic 
time variation of a are discussed in § V. 

3. Small m^. If <^ H , <f>, and therefore a, will be frozen at its present value. In this case, variation in a would 
only be possible if the value of <fi frozen in at different regions of space were different. Such models are discussed 
in § VI. 

While not every model will contain a minimum or be well modeled using a harmonic expansion about the minimum, 
it is straightforward to adapt the constraints on these three generic cases to a particular model. 
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III. RUNNING COUPLING CONSTANTS 



Before considering modifications to the QED Lagrangian in order to explain variation of the fine structure constant 
of the sort reported in Webb et al., we should consider whether the already known variation in coupling constants 
from predictions of quantum field theory might be responsible for this sort of deviation. Using the fine structure 
constant as an example, renormalization tells us that at high energies, the value of a must be different than at low 
energies since there are loop calculations involved. For one-loop diagrams, this "running" of a can be approximated 
as (e.g., [11]), for large momentum transfer Q 3> m e , 

£_ 

a ^ = /oN ' (9) 

1 _ -l^-ln ( — 



where m e is the mass of the electron and we have taken h = c = k = 1, or natural units. This running of a has indeed 
been observed experimentally, with Burkhardt and Pietrzyk [12] reporting 

- (91.1884 GeV) w 128.89 (10) 
a 

from measurements at Q corresponding to the mass of the Z boson. 
For Q « m e , the running of a can be approximated by [13] 

^~(^V. an 

a \m e J 

The resulting variation in a depends on the energy scales involved. To consider whether the variation would 
show up in current measurements of a, we consider the Si iv absorption line doublets used in [14], the most precise 
measurements that do not involve the potential difficulties of the many multiplet method. Si iv absorption lines come 
from the warm ionized interstellar medium (WIM), where the most energetic processes take place typically around 
10000 K [15]. 

Taking Q ~ 10000 K w 0.86 eV 

— ~3xl0" 12 . (12) 
a 

In addition to changes in a from renormalization group effects leading to a running coupling constant, there will 
also be finite-temperature corrections to any process involving a loop diagram. As a result, the Lamb shift will be 
altered and therefore the atomic spectrum of silicon. Depending on the temperature and on the atomic number of the 
atom involved, there are two corrections to atomic spectra that must be considered. Walsh [16] considers a correction 
to the Lamb shift dominant at high temperatures and at low atomic number, concluding that 

A£^ 300 — - (13) 

Z^(m e a) 6 

Using silicon [Z = 14) as an example, this correction would be 

AE w 2.6 x 10~ 22 eV. (14) 

Barton [17] considers a second correction, dominant at lower temperatures and higher atomic number over that 
calculated by Walsh, of approximately 



or again using the example of silicon, 

AE w 2.1 x 10~ 14 eV. (16) 

In both cases note that we have assumed a temperature of 10000 K. Since only the most energetic processes in the 
WIM take place on this energy scale, this is an upper bound on AE but might be a significant overestimate. 

We can also consider whether these effects could explain the many multiplet method results described in [1]. These 
also involve measurements of ions such as Mg II, Fe II, Al II, Al III, Zn II, and Cr II that occur in the WIM. Using the 
iron (Z = 24) line as an example, the correction to the Lamb shift is now 

AE w 4.8 x 10~ 14 eV (17) 

and the correction directly to a is again a few parts in 10 12 . Since the typical fine structure transitions in iron are of 
order 0.4 eV, these corrections are also far too small to cause the variation claimed in [1]. 
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IV. OSCILLATION 



As in § II, the equation of motion (8) allows for three typical behaviors of a massive scalar field in an expanding 
universe. For 3> H, the scalar field will oscillate with an amplitude that is damped by the term 2>H{t)<j). As 
an example, approximate the universe between a redshift of z = 4 and the present as matter-dominated, so that 
H(t) = This is a good approximation for 0.5 < z < 4 while as dark energy comes to dominate the Hubble 

parameter tends toward a constant value. Solving for the matter-dominated case, we find the amplitude drops off as 
Since the scale factor a(t) cx t 2 / 3 , the amplitude of oscillation in the recent past increases with higher redshift as 

A(z) cx (1 + z)^ (18) 

in the regime in which 3> H . 

Assume the value of a today is oscillating and that therefore the present value is at some random phase of the 
periodic cycle. The strongest current bound on a is [18] 

- = (0.2 ± 7.0) x 10" 16 yr _1 - (19) 
a 

With 95% confidence, a chosen at a random point from oscillation with amplitude A(z) and frequency / will be 

\a\ > 0.08A{z = 0)/. (20) 

Therefore, at the 2a level the atomic clock constraint is 

A(z = 0)/ < 9 x 10" 15 yr _1 a- (21) 

This bound does not apply to oscillations with T < 10~ 8 s <~ 10~ 2e /H , as atomic clocks cannot detect variation more 
rapid than the characteristic timescale. 

Further, one also must explain why measurements of the fine structure constant using isotopic abundances from 
the sub-nuclear reactor at the Oklo mine [3,5,19,20] yield a value so similar to the present value of a. The isotopic 
abundances found at the Oklo mine are consistent with a constant value of a that is sufficiently close to the present 
value as to reproduce the approximately 0.1 cV resonance for the interaction g| 9 Sm + n — >g2° Sm + 7. It has been 
shown that if a can be treated as constant during the lifetime of the reactor, the fine structure constant was within 
1.44 x 10~ 8 a of the present value while the reactor was operating [5]. 

For an oscillating model of a with sufficiently low frequency that the value of a can be treated as constant during 
the approximately (2.3 ±0.7) x 10 5 yr period of activity, the Oklo reactor could only operate at a time when the value 
of a were fixed so that the 0.1 eV resonance was reproduced. Assuming that the frequency is sufficiently high that 
there have been a large number of periods in the intervening 1.8 billion years, a slight change in the frequency causes 
the present value of a to sweep through multiple periods, so that we should again expect to be at a random phase 
today. As before, we require that at least 5% of possible phases would produce a present value within the tolerances 
set by the Oklo measurements. This is equivalent to a requirement at the 2a level that 

A(z = 0) < 1.4 x 10~ 7 a. (22) 

No work has been done that would rule out the possibility that in addition to being a good fit with a constant value 
of a, the isotopic abundances at Oklo might also be a good fit for a model in which the value of a and therefore the 
energy of the neutron capture would vary in a specified manner during the reactor's operation. 

Therefore, for oscillation at an unknown frequency the results from the Oklo mine strongly suggest a limitation 
on the amplitude of oscillation. We can further combine the atomic clock bound with observational constraints 
on frequency to produce an additional bound on the amplitude. For example, Varshalovich et al [21] measures a 
from the spectra of 1500 absorption systems between z = 2 and z = 4 (or a time window At « 0.15/ffo), f° r an 
average sampling frequency of f s <~ 10000ifo- If there were oscillation within that time window with amplitude 
A(z = 3) > 4 x 10~ 5 a, the oscillation would be detected in these results unless the frequency were greater than the 
Nyquist frequency, requiring 

T > 0.15/ffo or T < l/5000iJ - (23) 

The low frequency case is really monotonic time variation, and will be considered in § V. Combining these limits 
with (21), the allowable frequency range requires A < 3 x 10~ 8 ai, so we can conclude from a second set of constraints 
independent of Oklo that A(z = 3) < 4 x 10~ 5 a, or equivalently 
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i(z = 0)<5x lQ- 6 a. 



(24) 



As another example, we can consider whether oscillation could be consistent with the results of Webb et al. Clearly 
oscillation cannot be consistent simultaneously with both the variation claimed in Webb et al. and with the null 
result reported by Chand et al., but we will consider just the Webb results in isolation as an example of the difficulty 
of constructing an oscillatory model given present bounds. In that case, the amplitude of the oscillation today must 
be at least 

A(z = 0) w (0.81 x l(T 5 )a. (25) 

This is larger than the bound from (22) and the coincidence between the values at Oklo and today has a probability 
of 0.08%. The second (less likely) branch of Oklo would increase the probability to 1.6%, still requiring greater than 
a 2cr deviation from expectation. 

A(z = 0) is also larger than 5 x 10~ 6 a, and therefore the Varshalovich bounds on the frequency will apply. 
It is possible that future measurements will be sufficiently sensitive to evade the Varshalovich bounds on allowed 
frequencies. In that case, one can analyze these new measurements in a similar manner to the analysis performed on 
the Varshalovich data in order to produce applicable bounds. 

Returning to the example of Webb et al., we can now calculate the probability that an oscillating model with the 
given A(z = 0) is consistent with atomic clock measurements. Including the Varshalovich bounds on the frequency 
and assuming that a is at a random phase today, atomic clock limits on a today will be satisfied with probability 
1.1 x 10~ 4 . For the claimed variation to be the result of an oscillating model, both this bound and the apparently 
independent bound from Oklo must be satisfied, for a combined probability of 9 x 10~ 8 . 

In general, we can see that the constraints on models in which the fine structure is oscillating are very strong in 
the absence of additional physics that would both finely tune atomic clock measurements today and explain why the 
sub-nuclear reactor at Oklo was active during a period when the value of the fine structure constant was the same as 
it is today. At this point, an improvement of the precision with which A can be constrained obscrvationally would 
strengthen bounds on oscillating models. However, additional observations of absorption systems at the present level 
of precision will only serve as additional confirmation that oscillation does not exist in a region already ruled out by 
current observational and experimental bounds. 



V. MONOTONIC TIME VARIATION FROM SCALAR FIELDS 



In order to produce an a that is monotonically varying we can introduce a scalar field <j) and choose a function ot{<j)) 
such that the QED Lagrangian is (writing only the electron terms) 

Cqed = {d^f - V(<f>) + j e (i 7 „0" + y/Anaifl-fA,, - m e ) Ve - ^ (0mM„) - d u (eA^) f . (26) 

The equation of motion (6) requires VV" rj 3i?o , since small values of V" result in a frozen-in field 4> and large values 
result in oscillation with frequency ^> 3-ffo- 

Any function a{<p) is possible. In general, we can take a Taylor series 

a = a{<P) = a + Ai(0 - <fo) + ^A 2 {4> - 4> f + ■ ■ • , (27) 

where A„ = d n a/d<f) n and where we will take 0o = 0. How could one test for such a 0? Mass renormalization leads 
to a contribution to the masses of the proton and neutron that depends on a [22] : 

Aa 

Am„ = B„Aa w (0.63 MeV) (28) 

a 

Am n = B n Aa w (-0.13 MeV) — , (29) 

a 

defining B p and B n coupling constants for the proton and neutron, respectively. 

Dvali and Zaldarriaga [23] (see also Damour and Polyakov [8]) consider the example a((f>) = A</>, which leads to a 
Yukawa interaction 

£ = ■■■+ [BnXi^NN] . (30) 

N—p,n 
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FIG. 1. Feynman Diagram for the Yukawa interaction between two nucleons N and N' mediated by one <f> particle. 



They show this is a long-range force proportional to \\ (see Figure 1) but different for protons and neutrons, and 
that Eotvos experiments require 

< \- < 31 » 

We can similarly consider the coupling (A2 / A)(/> 2 N 'N . As shown in Figure 2, \2<f) 2 NN and by analogy the other 
higher-order terms in <j> also produce long-range forces, mediated by n (f> particles rather than just one. The A is the 
cutoff energy from renormalization, and has been brought outside of A2 in order to keep all of the A„ dimcnsionless. 




FIG. 2. Feynman Diagram for the cj> 2 interaction of two nucleons N and N' . 

In comparing the two long-range forces, the matrix element corresponding to the interaction mediated by n 
particles can be approximated in terms of the matrix element for the single <fi interaction as 

/ x \ 2 



Mn ~ A [\r) Ml ' (32) 

since the matrix element is proportional to the coupling coefficient squared. A is a constant with units of energy™ -1 , 
and can be approximated as using the center of mass energy y/s as y/s" 1 . The cutoff energy A is typically taken to 
be on the TeV scale. For a long-range interaction between nucleons separated by a distance r, y/s typically scales as 
1/r, and therefore the long-range interaction between two nucleons will have y/s <C A. The <f) n term produces a long 
range force with magnitude 

< 33 » 

/ s 2n-2 

The long range force is thus depressed by I ^ J . Regardless of whether the force is mediated by one 4> particle 

or many (f> particles, the contributions from each force will simply add since the external lines on each diagram arc 
identical. Therefore, we can conclude from Eotvos experiments that 

x x fVa\ V (da 1 r 2 a V 10" 10 

This is effectively an upper bound on the coupling of a to <j>, although its utility depends on the model, as the 
higher-order derivatives of a{4>) are only very weakly bounded. We will calculate some examples shortly. 

First, we must bound A(f> in order to use bounds on ^§ to bound Aa. For a flat universe, the kinetic energy 
density of our (f> field cannot be greater than the critical density. Dvali and Zaldarriaga [23] use this to show that 
Ai > ^ p . However, this is only a bound on the present value. Not only must the kinetic energy density fit this 
bound today, but it also must fit similar bounds in the distant past. In an expanding universe, the energy density of 
a type of energy with equation of state parameter 



G 



changes as a function of scale factor a as 



w =i = pm (SB) 



poza-^ 1+w \ (36) 



so that kinetic energy (w = 1) drops off as a~ 6 . Therefore, if there is a little kinetic energy now, we should expect a 
lot of kinetic energy when the universe was smaller, or equivalently at large redshift (a = 1/(1 + z)). For example, 
at a redshift of z = 1, the ratio of kinetic energy density to matter density would already be eight times larger than 
today; at a redshift of nine it will be 1000 times larger. The effects of even a small amount of kinetic energy today 
would be enormous in the early universe, as at high redshift it would dominate the universe. Among the effects we 
would have noticed would be a severe suppression of density fluctuations and a very large Aa. 

This does not prove that 4> cannot be rapidly changing today, but rather that, if <j) is changing, in addition to the 
kinetic energy of the field there must be a sufficiently large potential that < 0. Much of the "missing" energy in 
the universe, or the so-called dark energy, has w < 0, so it might be natural to postulate that the scalar quintessence 
field postulated to solve the dark energy problem is also responsible for variation in a. 

Having bounded A(j>, we can now return to our bounds on 0^ and use them to bound Aa. In order to highlight 
the limitations of the Eotvos constraints, let us choose as a toy model an a((f>) dominated by the <j> n term. We choose 

a = a Q +K(j) n , (37) 

with <p — at the present value of a. The bound from (34) thus yields 

VI d n a 10- 5 

n! 0^ M P ■ W 

In this case, 



so (38) becomes 



Since A(f> can be as large as Mp 3> \fs, the resulting bound of 



A " < (40, 



n-1 



does not provide a useful constraint for n > 1. 

Dvali and Zaldarriaga consider the specific case n = 1, for which ^ may be as large as I0~ 3 between the present 
value and observations at a redshift of 2. However the bound from the Oklo reactor is much more restrictive. For 
their model with constant <p and thus constant a, at the time of Oklo we should measure 

^) ~ 10- 6 . (42) 
a / Oklo 

Both branches of Oklo bounds rule out variation of this order. Indeed, at the present time the Oklo bounds are the 
strictest bounds on monotonic time variation. 

Bounds on violation of the equivalence principle have improved by a few orders of magnitude over the past decade 
or two [24]. The Satellite Test of the Equivalence Principle (STEP) has been proposed to improve the sensitivity of 
Eotvos experiments by five or six orders of magnitude [25] and therefore the bounds on da/ dtp may in the near future 
become more restrictive on a linear coupling than the measurements from Oklo. Particle accelerator results such as 
those discussed in § VI will be orders of magnitude weaker even than the Eotvos bounds. 

One might be tempted to conclude that because Eotvos experiments are only capable of bounding models with 
a = a + K(f> n for n = 1, models with n > 1 are capable of producing large variation in a between z = and 
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z = 1. Indeed, the best direct bounds on models with this type of time variation come from terrestrial bounds on and 
observations of the fine structure constant itself rather than from experimental bounds on associated effects such as 
violation of the equivalence principle. 

However, in order to produce such a model, the linear term in a(<p) must be very small and only the higher-order 
(^-derivatives of a contribute. In the absence of a symmetry (for example, a function of leading order </> 2 would require 
the symmetry <f> — > —(/>), this is an unnatural model. Even with a symmetry, the linear term will only disappear 
in an expansion about the minimum. Since the Eotvos experiments that constrain da/d(j) take place today, this 
requires that the value of a be very near a minimum at present. Since we are dealing with monotonic oscillation, the 
requirement is for the value of a to be very close to the minimum today but to be far from the minimum both in the 
recent past and in the near future. Therefore just like the oscillating models discussed in § IV, such a model requires 
a very improbable coincidence. 

We do know that something fundamental has changed about the universe between a redshift of 3.5 and today. The 
universe has gone from decelerating to accelerating, and has gone from matter-dominated to a universe dominated 
by a dark energy with negative pressure. Therefore, one might naturally consider a model in which either a or <f> is a 
function of the composition of the universe, perhaps a non-minimally coupled scalar field (e.g., the Brans-Dicke [26] 
model with a scalar field coupling to the Ricci scalar). 

The coincidence between the timing of dark energy and of the reported variation in a from Webb et al. suggests 
that the two might be related. Wetterich [27] suggests that such a model might satisfy both variation on the order 
reported in Webb et al. and the Oklo bounds. In order to produce the large variation between z = 0.1 and z = 2 
compared to the variation between today and z = 0.1, the dependence of a on dark energy must be similar to a step 
function. 

Assume that the dark energy is a time- varying quintessence field Q rather than a cosmological constant. Then the 
kind of dependence one is looking for might be of the form 



Letting, as an example, Qq = 0.1, a would take on a constant value in the distant past as well as a constant value 
today, but the constant value would be different in those two epochs. Similar ideas have been discussed by Barrow 
[28] among others. However, no such model that fits both the variation reported by Webb et al. and all of the other 
measured bounds, from not only Oklo but also the CMB [29] and nucleosynthesis [30], has yet been worked out in 
detail. While Barrow [31] suggests that a model with some features in common with the one proposed here will have 
the correct asymptotic features, it is not at all clear that the boundary conditions can be matched in such a way as 
to actually fit current measured limits. 

To this point, a model that both is consistent with the claimed variation in [1] and arises naturally from current 
theory has not been forthcoming. Additionally, for time variation the recent results reported by Chand et al. are in 
direct conflict with those of Webb et al. Any model in which the variation of a is time requires that independent 
measurements of a over the same range of redshift must produce the same result. Therefore an additional requirement 
of any model that produces variation of the sort claimed by Webb et al. by introducing some a(Q) is that either the 
Chand et al. results are shown to be erroneous or Qq is shown to have already taken on its present local value in 
some regions of the universe at a redshift of 2.3, which would be inconsistent with inflation. 



Another logical possibility for producing a different value of a at z — 1 than at z = is that rather than changing 
over time, a undergoes spatial variation. This is perhaps the simplest way of removing the Oklo constraint, as spatial 
variation could produce a different value of a at large distance (or, equivalently, at large red shift) while retaining 
the same constant value on Earth in the distant past. Oscillating models will run into similar constraints as § IV, 
although only the bound (21) will be pertinent. 

In order for spatial variation to occur, there must be (at least) two values of a for which it will remain constant. 
One possibility, which we will not consider here, is that a light field, such as that postulated by "tracker" quintessence 
models, is stuck as in [32]. For small V and V" <C H 2 , <j> will be small. Such models typically produce small V and 
V" by requiring the field to couple to almost nothing else; since a must be a function of our scalar field it will couple 
to almost everything. 

The other possibility is a relatively heavy field in which the potential V(a) has at least two local minima, cor- 
responding to the values of a locally and at large distances. In the early universe, the temperature is very high, 
and for total energy much greater than the potential, <p can take on any value. As the universe expands and cools, 
at some point (f> falls into one of these local minima and cannot get out of the potential well. However, the local 



-n, 



ii. 




(43) 



VI. SPATIAL VARIATION 
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value determines which well <f> will fall into, and thus different regions of the universe can randomly fall into different 
potential wells. Therefore, it is possible that we are in one potential well with one value of a and much of the universe 
is in another potential well with a different value. 

As in the previous section, we can equivalently consider a real scalar field <ft such that a is a function of <p. We 
select the minima to be ±<r. We can then write the Lagrangian as 

C = \(d^) 2 -\\(tf-^)\ (44) 

where A is typically of order unity. 

At high temperatures, d^cp 3> § 2 ~ 0-2 , and therefore the expectation value (<f>) = 0. However, at lower temperatures 
<f> will take on a vacuum expectation value of either (</>) = a or (cf>) = — cr, spontaneously breaking the 2,2 symmetry 
under <f> — ► — <fi. We might imagine that the local vacuum has ((f)) — a and that the vacuum containing absorption 
systems measured by Murphy et al. has ((f)) = — a. Since (f>(x^) must be continuous, even at very low energies there 
will be a region of space called a domain wall for which = 0. The major problem that we shall run into will be the 
total energy density of these walls. 

A scalar field has stress tensor 



Tfxv = d^d v <t> - g^C. (45) 

Because the domain wall is in a region of false vacuum, instead of = the potential term of the Lagrangian yields 
[33] 

T \ = lAa 4 cosh- 4 (A 1 / 2 ^) (46) 

for a wall located in the x — y plane centered at z = 0. 

So, what is the energy of the domain wall network? For small values of z, i.e. values at and around the domain 
wall, 

T° « (47) 

As a lower bound on the wall energy density p w , consider the case in which there is only one domain wall in the 
observable universe. Then, p w will be the surface energy density multiplied by the wall thickness and divided by the 
horizon size r, or 

Pw ~ ^\<7 4 z w /r. (48) 



From (46), 

1 



Note that we have set U = c = 1. For a flat universe (i.e. k=0) such as ours, 



(49) 



wrm (50) 



Since we can expect this symmetry to have been broken during the radiation-dominated era (more on this shortly), 
the energy density of the universe is p w T 4 , and by the Friedmann equation 



Further, the symmetry is broken when the thermal energy drops low enough that <f> gets trapped in one of the minima 
of the double well potential in (44), or when T w crA 1//4 . 
Therefore, the total energy density of the domain wall is 

1 Act 4 /8tt a 1 



2A 1 /2 (T V 3 M P 
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Pw « (52) 

Since the universe is flat, the critical density at that time is very nearly the radiation density 

p c wp r «T 4 «A<t 4 . 

So, at the time of spontaneous symmetry breaking (SSB), the ratio of the energy density in domain walls to the 
critical density will be 

^ssb « A l/2 Mp ■ ( 53 ) 

Since a is an electroweak coupling constant, this symmetry cannot be broken until at or after electroweak decoupling 
at approximately 300 GeV [33]. Since Mp w 10 19 GeV, for A 3> 10~ 31 (which is shown to be the case below), il w <C 1 
and therefore we can treat domain walls as negligible in the Fricdmann equation describing the expanding universe. 
WMAP bounds using the Integrated Sachs- Wolfe effect confirm that fl w < 10~ 4 today [34] compared to TO = 0.28, 
so domain walls can be ignored in the Friedmann equation for the entire history of the universe. 

The energy density of a domain wall network is known to scale with the size of the universe as approximately 
p w cx a -1,44 [35]. This will again underestimate p w , since a solitary wall has p w cx a -1 and the additional dropoff is 
due to coarsening, based on interactions between different causally connected walls. 

Let the scale factor at the time of symmetry breaking be a w = 1/(1 + z w ), where today a = 1. Then, 

n W0 ~n WSSB p Esssb^ (54) 

\ a wJ Pc 

During the radiation-dominated era, i.e. z 3> 3500, p c cx a -4 and during the matter-dominated era, z <C 3500, 
p c cx a~ 3 . (The effects of A-domination since a redshift of 0.5 or so are negligible). So, for 1 + z w > 3500, 



If instead the symmetry is broken during the matter dominated era, 

Qw ° K AV^ (1 + ^ )1 ' 56 ' (56) 

so (55) is a lower bound on fl Wo . 

In general, the scale factor is inversely proportional to the temperature, a cx T _1 . Wc know that the symmetry is 
broken at approximately a temperature of a. Therefore, we can conclude that 

l + « ro w=-^— . (57) 

J CMB 

where Tcmb is the current temperature of the universe, 0.000235 eV. Substituting, we conclude that 

CT 3.56 

Qwo ~ (5.65 x 10-7)AV2Mp' (58) 

with a in eV. 

Q w < 10 -4 today, so substituting we conclude that 

(jA -0.140 < 103 ke y ( 5Q ) 

Is this upper bound valid in all cases? If the formation of walls occurs prior to inflation, the walls can be inflated 
far away, similar to the solution of the monopolc problem. It might then simply be an accident that one happens to 
lie in our observable universe. However, as mentioned above, the electroweak decoupling only occurs around 300 GeV, 
well after inflation. Further, if the domain wall were inflated out of our observable universe, our observable universe 
would all be in the same domain. Therefore, this model could not produce variation in the fine structure constant. 

Another way to prevent domain wall problems is restore the symmetry at a lower temperature. At first, this might 
also seem problematic, since we want the spatial variation laid down initially to remain today. Indeed, such a model 
does not properly fit under the spatial variation heading at the start of this section. Due to causality, we cannot 
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observe the present value of a at large distances, and therefore we cannot rule out the possibility that the value of a 
is different at present in gas clouds large distances than it was when light we observe today was passing through those 
regions. However, since the observations of Webb et al. include systems at z = 0.5, WMAP bounds on Cl w using the 
Integrated Sachs- Wolfe effect since z = 2 will continue to constrain such models. 

A third option would be for there to be a slight asymmetry, such that the two local minima are unequal. In that 
case, one vacuum would be preferred with regions of higher-density vacuum shrinking and eventually becoming black 
holes. However, this effect only saves the model from becoming wall-dominated when these regions start to disappear, 
so a model that postulates two regions of different vacuum today or in the recent past is incompatible. 

By introducing a scalar field </>, we have also introduced a particle <p. What is the mass of this particle? As in (44), 

C = \(d^) 2 -\\(tf-a 2 ) 2 . (60) 



The mass at the present value (j) = a will be = y/V"{a), or 

m = A5oV2. (61) 

(As an aside, notice that had we allowed a complex scalar field it would also yield a particle of mass Aa\/2 along with 
a Goldstone boson.) 

We can constrain \o\f2 for a model of the sort that might fit measurable variation by requiring that the symmetry 
is broken after electroweak decoupling and before z — 3, or 

300 GcV > T w crA 1/4 > 10 K. (62) 

Combining this with the WMAP bound on domain wall density from (59), we find 

\?aV2 < 5.7 x 10~ 9 eV. (63) 

Because there is no lower bound on the wall density, there is no direct upper bound on A. As A gets large, constraints 
on crA 1 / 4 require that a approaches zero and the Lagrangian simply describes A</> 4 theory. Having fixed A on some 
energy scale M, the effective A will be [13] 

= 1- (3A/16^ 2 )log(p/M)' (64) 

When (3A/167T 2 ) log(p/M) = 1, the perturbation theory breaks down. We require that the theory be valid for 
p < 300 GeV, the range on which accelerator experiments may be able to test it. If we fix A at M = 91 GeV, the 
energy scale of the relevant L3 measurements discussed below, this requires A < 41.6, and therefore 

m < 760 GeV. (65) 

Because variation in a depends on the expectation value of 4> rather than on the presence of the associated particle, 
it is acceptable for the mass of the particle to be far above the background temperature at the time of spontaneous 
symmetry breaking. However, since in the near future there will be experimental results at energies greater than even 
the largest possible m^, let us consider the case of a particle with mass less than 91 GeV and ask whether it would 
have been observed. 

The QED Lagrangian includes a term 



-QED 



= eV4^7 AI A Al e, (66) 



where e is the electron field. Since we have postulated a cf> — > — <fi symmetry, we can expand a = ao + Since 
Aa = ai — ao is small compared to ao, we can approximate i/a y / oo + Over a change in <j) of la a changes 

by Aa. While we do not know the function a(4>), we can estimate the resulting bound by assuming a constant 
Then, 

da^ _ Mrf) = Aa 

d(j> 2 ~ A(ct) 2 ) 2a 2 ' { ' 

So, the coupling eA^e splits into two branches, 
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£=■■■ + V^e (l + \ (0 J 2 j y^e- (68) 

Known electron-positron decays begin with e + e~ decaying first to a virtual photon, and then to real products. We 
have now introduced a different decay in which e+e~ — > 70</>. No <j) particles have been detected, so this decay should 
appear to be a single-photon decay forbidden by conservation of momentum in the center of mass frame. However, the 
coupling producing particles is a five-point interaction with mass dimension six, and is therefore Planck suppressed. 
Therefore, even the best L3 bound [36] of 7 e +e-->7^/7other — ^-^ x 10 -6 cannot be used to constrain this theory. 

It is also possible that the Zi symmetry in the effective Lagrangian is explicitly broken instead of spontaneously 
broken, and that there is a coupling e<f>A^e. This would yield two unequal minima and the corresponding problems 
discussed earlier. We can repeat the calculation above using a coupling to a single <fi, but even so the mass dimension 
will be five and the four-point coupling will be Planck suppressed. 

So, this model is consistent with every observational and experimental bound aside from the many multiplet results. 
Webb et al. report 

.(0.2 < z < 3.7) - a(z = 0) = ( _ Q ^ ± Q1Q) x ^ (gg) 



while Chand et al. report 



a(z = 0) 



a(0.4 < z < 2.3) -a(z = 0) _ (wu; n n; .. . m , 



a(z = 0) 



= (-0.06 ±0.06) x 10"°. (70) 



The two measurements are inconsistent with one another if a is spatially uniform, but closer examination is required 
if a is spatially varying. In particular, one might examine whether the discrepancy can be explained by the fact 
that the two groups examined different regions of the sky. The Chand et al. results are based on a sample from the 
Southern Hemisphere only, whereas Webb et al. consider systems in both the Southern and Northern Hemisphere. For 
simplicity, we divide the Webb et al. data into Southern only and Northern only. (A more sophisticated study might 
break the Webb et al. dataset along a different boundary, but our choice will suffice for the purposes of illustration.) 
For 96 quasars in the Northern sample, 



a(0.2 < z < 3.7) - a(z = 0) 

"0 = 0) /North 

while for 32 quasars in the Southern sample, 

'a(0.2 < z < 3.7) - a(z = 0) s 



a(z = 0) /South 



(-0.66 ± 0.12) x 10"°, (71) 



(-0.36 ± 0.19) x 10" 5 . (72) 



Therefore, this model allows for the possibility that both Chand et al. and Webb et al. are consistent despite 
appearance to the contrary. Additional observations by both groups in both Hemispheres would allow a determination 
of whether the two sets of analysis agree on the value of a when looking at an identical patch of sky. If so, and if 
both groups agree that there is a different value in parts of the Northern and Southern hemispheres, this would be 
strong evidence for spatial variation. Such a test could be accomplished with an additional sample of approximately 
50 quasars and would be capable of either ruling out or confirming whether such a model can explain the discrepancy. 
Alternatively, this test might instead show a systematic difference between the results of the two groups in identical 
patches of sky, which would suggest that the apparent discrepancy should be resolved through a re-examination of 
systematic errors rather than through a new theoretical model. 



VII. DISCUSSION 



Observations from distant absorption systems suggest that the fine structure constant may not be constant, but 
rather varying as a function of distance or time. While other measurements suggest the value of a might instead 
be constant, in this paper we have considered features generic to models that predict spatial or time variation. In 
particular, we have considered models in which a varies as a function of some scalar field, since such models could be 
motivated either by string theory or more recently by cosmological attempts to solve the dark energy problem. Scalar 
fields can cause any of three classes of variation of a, considered in § II. 
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• Oscillating a. If a is oscillating rapidly in the recent past, the characteristic change in a is constrained by 
measurements of a from atomic clocks today. For oscillation of amplitude A and frequency /, this requires 



A(z = 0)/<9x l(T 15 a 



(73) 



except at frequencies greater than 100 MHz. For sufficiently large amplitude, it would also be a coincidence 
that the value of a at Oklo is very nearly the same as a today. This coincidence requires 



• Monotonic time variation of a. If a is changing monotonically in the recent past, the strongest bound 
comes from the Oklo sub-nuclear reactor. For the simplest model with a{4>) = a a + K<p and with constant <p 
than the Oklo bound 



is stronger than any phcnomcnological bounds. Projected Eotvos-type experiments such as STEP may improve 
this bound. 

Eotvos bounds on higher-order couplings of a to do not provide useful constraints. 

• Spatial variation of a. If a is changing with distance rather than with time, there are two possibilities. A 
light, minimally coupled field could result in a monotonic a{x l ), but it is difficult to have a minimally coupled 
field that also causes a change in the ubiquitous fine structure constant. Or, a heavier field with a discrete 
symmetry can undergo spontaneous symmetry breaking, resulting in different values of a in different regions of 
space. 

We have shown that the Chand et al. and Webb et al. results agree in the Southern Hemisphere, where their 
sampled regions overlap. Both are consistent with no variation in a from the terrestrial value. However, Webb et al. 
also include a sample from the Northern Hemisphere which does show a statistically significant deviation from the 
terrestrial value. If a is uniform, then the two results are inconsistent. However, we have found that a model in which 
a has different values in different domains (with one value for the Earth and Southern Hemisphere and another for 
the Northern Hemisphere) is consistent with both measurements and with all other current constraints. Alternatively, 
the observed difference between Southern and Northern Hemisphere may be due to a systematic bias that has not yet 
been identified. Further observations should clear up the situation. 

Another possible explanation for spatial variation of coupling constants is the recent "Chameleon" model in which a 
scalar field can have a mass dependent upon the local matter density [37]. Such a model would allow the 4> field to take 
on one value terrestrially and another value in gas clouds, and therefore allows variation that cannot be constrained 
by Oklo or atomic clocks. However, unlike the spatial variation case discussed in § VI, there is no artifact such as a 
domain wall that would restrict such a model. The only additional constraint that may prove useful would be from 
nucleosynthesis, but there are no restrictions on the form of a (</>), and therefore it should be possible to satisfy this 
constraint. Similar behavior also may be obtained without a direct coupling to the density, as proposed by Barrow 
and Mota [31,38]. To first approximation, gas clouds in different regions of the sky should have the same density, 
and therefore we should not observe spatial variation of a between gas clouds in different patches of sky. However, 
since the desired effect on a is quite small, it maybe that only a tiny but systematic difference between the densities 
of gas clouds in one region of sky and in another is required to produce differences in a of a few parts per million. 
Experiments scheduled for the next decade or two should be sufficiently sensitive to test predicted variation of the 
gravitational constant G according to this model [37] . 

Certainly if spatial variation of a were confirmed in the near future, this would provide an immediate and strong 
constraint on allowable models for unified theories. However, it should be stressed that any scalar field introduced 
that would cause the variation of a is subject to the bounds discussed in this paper. Therefore, improving current 
constraints on the variation of the fine structure constant and other coupling constants will continue to restrict 
particular models for unification. As there are at present very few measurements of observations with the potential 
to restrict such theories, it is even that much more critical to resolve the issues with the two many multiplet results 
and to continue to improve these bounds. 



A(z = 0) < 1.4 x lO^a. 



(74) 



|d| < 10 _17 a yr 



,-i 



(75) 
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